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ABSTRACT 

We present an investigation of the scale-dependence of bias described by the 
linear model: {6p{x.)/p)g = b{6p(x)/p)m, b being the bias parameter, and p{x)g and 
p(x)m are the galaxy number density and mass density, respectively. Using a discrete 
wavelet decomposition, we show that the behavior of bias scale-dependence cannot be 
described by one parameter b. In the linear bias model the scale-dependence should be 

~(n) 

measured by the j-spectra of wavelet-coefficient-represented bias parameters bj and 

(n) ~{n) 

bj , n being positive integers. Because bj with different n are independent from each 
other, a systematic analysis of the j-spectra of 6^"^ and bj^^ is necessary. 

We performed a j-spectrum analysis for samples of elliptical and lenticular 
(EL), and spiral (SP) galaxies listed in the APM bright galaxy catalog. We found 
that, for statistics of two-point correlation functions or DWT power spectrum, the 
scale-independence holds within 1 a. However, the bias scale-dependence becomes 
substantial when phase-sensitive statistics (e.g. 5^"^ with n > 2 or b^^^) are applied. 
These results indicate that the bias scale-dependence has the same origin as the 
non-Gaussianity of galaxy distributions. This is generally consistent with the 
explanation that the bias scale-dependence originated from non-linear and non-local 
relationship between galaxy formation and their environment. 
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1. Introduction 

Bias is introduced to reconcile the amplitude of fluctuations inferred from clustering of 
galaxies with that derived from mass distributions. Bright galaxies seem to have a stronger 
clustering than that of the underlying mass. Therefore, it is generally believed that galaxies are 
biased tracers of the mass density field, i.e. it follows the clustering of (dark) matter but with an 
enhanced amplitude. 

Generally, bias is phenomenologically modeled by a linear relation as 

5(x)g = b6{^)m. (1) 

where 5{'x)g = [n(x) — n]/n, 5{x)m = [/o(x) — p]/p, /3(x) and n(x) are, respectively, the galaxy 
number density distribution and density field of dark matter, and p and n being the average of 
p(x) and n(x). The phenomenological bias parameter b is assumed to be constant, but may be 
different for different galaxy types, say, beariy, hate for early and late types of galaxies. 

Eq.(l) and its variants are widely used in the determination of cosmological parameters from 
samples of redshift surveys of galaxies. However, we have really very little idea about which the 
bias parameter b should be. In fact, both bias (6 > 1) and anti-bias (6 < 1) are employed in 
current data analysis (e.g. Mo, Jing & White 1996.) This prevents unambiguous measures of 
cosmological parameters, giving only bias-contaminated results. 

Theoretically, the physical mechanism responsible for relation (1) is far from clear. The 

first analytic model of bias, in which objects arc identified with high peaks or collapsed halos of 
the density field, offers a plausible explanation of the bias of galaxy clusters - the correlation 
amplitudes of clusters are strong functions of cluster richness (Kaiser 1984). In this case, bias is 
mainly caused by the mass of collapsed halos. The larger the mass of the halo, the higher the 
richness of the cluster. 

Yet, the formation of galaxies doesn't depend only on the mass, or local mass density, of 
collapsed halos, but is substantially modulated by various environmental effects, such as the 
suppression of star formation in neighboring protogalaxies (Rees 1985), stimulating the formation 
of nearby galaxies (Dekel & Rees 1987), dynamical friction effects (Couchman & Carlberg 

1992) , etc. All these environmental effects are beyond local density, and lead to a non-local 
relation between the number density of galaxies and the background mass field (Bower et al. 

1993) . Moreover, the rate of star formation is most likely non- linearly dependent on local mass 
density. A common result of the non-local and/or non-linear relation between 5g and 5^ is the 
scale-dependence of parameter b. Therefore, to have a deep understanding of the mechanism of 
galaxy bias, searching for the scale-dependence of parameter b is necessary (Coles 1993, Catelan 
et al. 1994.) 

So far, the results of detecting b scale-dependence are quite scattered. For instance, the values 
of the bias-contaminated density parameter (3 = Cl'^-^/b are found to be in the range of 0.4 - 1, and 
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equal to about ~ 0.5 — 0.6 at Gaussian smoothing scales of 3-6 h~^Mpc, and ~ 1 on scales of ~ 12 
h~^Mpc (e.g. Dckcl, Burstein & White 1996). This is, b is probably scale-dependent from 6 to 12 
h^^ Mpc. On the other hand, some studies conclude that for galaxies of all types and luminosities 
the scale dependence of the bias parameter is weak (e.g. Kauffmann, Nusser & Steinmetz 1997). 
Why different detections gave different conclusions? This question motivated us to study the 
physics included in Eq.(l). 

In the first part of this paper, we show that even in linear bias model, the behavior of 
scale-dependence cannot be described by one parameter b, but by a series of j-spectra of 
wavelet-coefficient-represented bias parameters b^^^ and bj^\ n being positive integers. The 

~(n) 

j-spectra of bj are all statistically independent. Different detections may actually measure 

~(n) 

different parameters bj . Therefore, it should not be surprised that some detections arc positive, 
and some negative. The different behavior of bias scale-dependence given by different detection 
does not cause confusion, but may greatly be helpful to reveal the physics behind the bias model 
(1). Therefore, to have a complete picture of bias scale-dependence a systematic analysis of the 
j-spectra of 6^"^ and 6^"^ is necessary. 

In the second part, we performed a systematic detection of the bias scale-dependence with 
the samples of galaxies listed in the APM bright galaxies catalog (APM-BGC). This analysis 
shows that for the APM-BGC sample, the scale-independence approximately holds if only the 
two-point correlation function and power spectrum are involved, while the scale-dependence 
becomes substantial when higher order or phase-sensitive statistics involved. This result is worth 
to constrain models of hierarchical clustering of galaxies. 

The paper is organized as follows. In §2 we briefly introduce a method of the space-scale 
decomposition based on discrete wavelet transform (DWT) analysis. With this method, various 
statistics of measuring the bias scale dependence are developed and presented. §3 describes the 
sample to be analyzed. In §4 we apply the DWT method to analyze the galaxy sample of the 
APM bright galaxies. The implications of these results are discussed in §5. 

2. Linear bias model in DWT representation 

2.1. The discrete wavelet transform 

Let us briefly introduce the DWT analysis of large scale structures, for the details referring to 
(Fang &; Pando 1997, Pando & Fang 1996, 1998.) We consider here a 1-D mass density distribution 
p{x) or contrast 5{x) = [p{x) — p]/p, which are mathematically random fields over a spatial range 
< X < L. It is not difficult to extend all results developed in this section into 2-D and 3-D 
because the DWT bases for higher dimension can be constructed by a direct product of 1-D bases. 
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Like the Fourier expansion of the field S{x), the DWT expansion of the field S{x) is given by 

CO 2^ -1 

Six)=Y,Y.^M^) (2) 
j=o 1=0 

where tpj^i(x), j = 0, 1, .... I = 0...2-' — 1 are the bases of the DWT. Because these bases are 
orthogonal and complete, the wavelet function coefficient (WFC), wj^i, is computed by 

= J 5{x)iljj^i{x)dx. (3) 

The wavelet transform bases ipj,i{x) are generated from the basic wavelet ip{x/L) by a dilation 
2^, and a translation I, i.e. 

V',,K^)=(^1 H^'x/L-l). (4) 

The basic wavelet tp is designed to be continuous, admissible and localized. Unlike the Fourier 
bases exjp{i2Tmx / L) , which are non-local in physical space, the wavelet bases localized 
in both physical space and Fourier (scale) space. In physical space, ^pj^i{x) is centered at position 
lL/2^, and in Fourier space, it is centered at wavenumber 27r x 2^/L. Therefore, the DWT 
decomposes the density fluctuating field 5{x) into domains j, I in phase space, and for each basis 
the corresponding area in the phase space is as small as that allowed by the uncertainty principle. 
WFC Wj^i and its intensity \wj^i\'^ describe, respectively, the fluctuation of density and its power 
on scale at position IL/2K 

In order to have reasonable statistics, the cosmic density field is usually assumed to be ergodic: 
the average over an ensemble is equal to the spatial average taken over one realization. This is the 
so-called "fair sample hypothesis" (Peebles 1980). A homogeneous Gaussian field with continuous 
spectrum is certainly ergodic (Adler 1981). In some non-Gaussian cases, such as homogeneous 
and isotropic turbulence (Vanmarke, 1983), ergodicity also approximately holds. Roughly, the 
ergodic hypothesis is reasonable if spatial correlations are decreasing sufficiently rapidly with 
increasing separation. The volumes separated with distances larger than the correlation length 
are approximately statistically independent. In this case, WFCs, Wj^i, at different I from one 
realization of 5{x) can be employed as statistically independent measurement. Thus, the WFCs 
of 5{x) at different I on a given j, iD^y, form an ensemble of the WFCs on the scale j. In other 
words, when the "fair sample hypothesis" holds, an average over ensemble can be fairly estimated 
by average over I, i.e. {wj^i) ~ {l/2^)Y^fj^^ Wj^u where (...) denotes ensemble average. 

2.2. Reconstruction of density field 

Using the completeness of the DWT basis, one can reconstruct the original density field from 
the coefficient Wj^i. To achieve this, DWT analysis employs another set of functions consisting 
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of the so-called scaling functions, (f)j^i, which are generated from the basic scaling (j){x/L) by a 
dilation 2-', and a translation Z, i.e. 

4>jA^)=[j] 0(2V^-O- (5) 

The basic scaling cj) is essentially a window function with width x/L = 1. Thus, the scaling 
functions 4>j,i{x) are also windows, but with width {1/2^)L, and centered at . The scaling 
functions <j)j,i{x) are orthogonal with respect to the index Z, but not for j. This is a common 
property of window functions, which can be orthogonal in physical space, but not in Fourier space. 

For Daubechies wavelets, the basic wavelet and the basic scaling are related by recursive 
equations as (Daubechies 1992) 

<t>{^/L) =Y.iai<P{2x/L-l) 

^{xjL) = Y.I bi4>{2x/L + 1) ^ ' 

where coefficients a/ and bi are different for different wavelet. In this paper, we use the Daubechies 
4 wavelet (D4), for which oq = (1 + v^)/4, ai = (3 + V3)/4, as = (3 - \/3)/4, 03 = (1 - V^)/i. 

From Eq.(6), one can show that the scaling functions (pj,i{x) are always orthogonal to the 
wavelet bases ^j'^i'{x) if j < /, i.e. 



(t>j,i{x)ipji^li{x)dx = 0, for j <j . (7) 

Therefore, (f>j,i{x) can be expressed by 'tpji^ii{x) as 

00 2^' -1 j-l 2J'-1 

^jA^) Cji.3'i''^j',i'{x) =^Y1 (^ji;j'i'i^j',i'ix). (8) 

j'=0 l'=0 j'=0 l'=0 

The coefficients Cji-j'ii = J (f)j^i{x)il>ji^i'{x)dx can be determined from ai and bi. 
Using (f)j^i{x), we construct a density field on scale j as 

2^-1 

1=0 

where Wj^i is called the scaling function coefficient (SFC) given by 

= f p{x)(j)j,iix)dx. (10) 
Jo 

Since the scaling function (i)j,i{x) is window- like, the coefficient Wj^i is actually a "count-in-cell" in 
a window on scale j at position I. 
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Using Eqs.(2), (8), (9) and (10), one can find 

j-l 2i'-l 

p'{x) = pj2 '^j',i''^3',i'i^) + p- (11) 

ji=0 11=0 

Namely, fP{x) contains all terms of density fluctuations Wji^ii'tpji^ii{x) of / < j, but not terms of 
f > j- Prom Eqs.(2) and (11), we have 

p{x) = pP{x) +pYY. ^j',i'i^j',i'{x). (12) 

j'=j l'=0 

One can also define the smoothed density contrasts on scale j to be 

s'{x) ^ ^^^-^ = E E ^j',i'i^j'A^)- (13) 

P j'=0 l'=0 

Eqs.(12) and (13) show that p^{x) and 5'^{x) are the density field smoothed on scale j. Using 
Eqs.(12) and (13), one can construct the density field />'(x) or 5^{x) on finer and finer scales 
by WFCs Wj^i till to the precision of the original field. Since the sets of bases Vj,; and (pj^i are 
complete, the original field can be reconstructed without lost information. 



2.3. Statistical description of scale dependence of bias 

We subject the DWT decomposition on Eq.(l) which gives 

(^j,i)ff = Kwj,i)m- (14) 

Statistically, Eq.(14) means that the one-point distributions of the WFCs for galaxies should 
be proportional to the matter, and the proportional coefficient h is j-independent. Thus, linear 
bias models Eq.(l) requires that various statistics based on the WFC one-point distributions for 
galaxies and matter should satisfy the linear relations given by Eq.(14). 

If we define the cumulant moment of the one-point distribution of Wj^i as 

^;=(-^^) = ^E-"^ (15) 

z=o 

the model (1) is actually to require that the ratios {CJ)g/{C'^)m = {{'Wj^i)g) / {{wj^i)J^) should be 
j(scale)-independent for all order n. Namely, 



r(n) ^ (C")g 



(16) 
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~(n) 

should be j-independent. In other words, the j-spectra of bj ' are flat. For non-Gaussian fields, 

~(n) 

the parameters bj ' for different n generally are statistically independent. Therefore, the behavior 

~(n\ 

of bias scale-dependence cannot be described by one parameter 6, but the j-spectra of 6^- ' for all 
n. 

Prom Eqs.(13) and (14), one has 

8^{x)g = b8\x)^. (17) 

Eq.(17) can be rewritten as 

f^{x)g -Pg = l)\p'(x)„^ - p„,] (18) 

where constant h' = hpg/p„i- Considering that / (j)j^i{x)dx is independent of index Z, the DWT 
decomposition of Eq.(18) gives 

{Wj^l - Wj)g = b'iWj^i - Wj)m (19) 

where Wj = (1/2-') Y^^q^ 'Wj,i is the mean of the SFCs on scale j. Similar to Eq.(16), one can 
design statistics fc^"^ as follows 



3-1 
=0 

(n) 



^(n) _ ^t=u ~j/y ^2{)^ 



If bias is scale-independent, the j-spectra of 6^ ' should also be flat. 

~(n) 

In a word, the scale-dependence of bias is described by the j spectra of parameters h - and 

(n) 

bj . It is possible that some j-spectra are flat, but others non-flat. Therefore, to study bias 
scale-dependence, it is necessary to take a systematic detection of the j-spectra of and b^J^^ . 



2.4. Non-Gaussianity and spectra of bj and bj 

A physical reason of choosing 6^"^ and 6^"^ to detect scale-dependence of bias is from 
non-Gaussianity of mass density fields. 

The structure formation due to gravitational clustering can roughly be divided into three 
stages: linear, quasilinear and fully developed nonlinear. It is obvious from eq.(l) that the linear 
evolution of S{x)g and 5{x)m do not cause bias scalc-dcpcndcncc. It is generally believed that 
the evolution of cosmic clustering on very large scales is probably still remaining in linear stage 
at present. Therefore, the bias parameter should be scale-independent on very large scales if the 
initial b is scale-independent. 

On scales from about 5-10 h~^ Mpc, the cosmic gravitational clustering may also be not 
fully developed yet, but partially in the quasilinear regime. In this case, the power spectrum of 
density perturbations is not significantly different from linear power spectrum, but is characterized 
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by power transfer via mode coupling. The power transfer among perturbations on different scales 
will give rise to a scale dependence of bias. On the other hand, power transfer also leads to 
non-Gaussianity of the density fields. Therefore, the bias scale-dependence in the range of 5 
- 10 h~i Mpc may not easily be detected by methods insensitive to non-Gaussianity, such as 
power spectrum or two-point correlation function, but by statistics sensitive to non-Gaussianity. 
The statistics and 6^"^ are non-Gaussian sensitive, and therefore, suitable to detect the bias 
scale-dependence caused by non-Gaussian process. 

It has been shown that the n = 2 statistics of (wj^i), i.e. (t| = C|, essentially is the Fourier 
power spectrum. Therefore a Gaussian random field can be completely described by its DWT 
power spectrum (t| = Cj (Pando &; Fang, 1998.) Higher order (n > 2) statistics of (wj^i)"' are 
sensitive to non-Gaussianity of the density fields. For instance, the n = 3, 4, statistics of {wj^iY 
can be rewritten as 

Therefore, Sj and Kj arc essentially the same as the statistics of skewness and kurtosis. This 
point can directly be seen from the definition of the skewness and kurtosis of density contrast 
distribution 5{x). They are 



S = A,C^ K = —,C^ - 3 (22) 



where C" = (1/L) /^f and = C^. Therefore, Sj and Kj are the j-spectra of S and K. 

For Gaussian fields ?){x), all the cumulant moments C" or are equal to zero, only except 
second order or (Fang h Pando 1997.) Thus, all statistics 6^"^ of n > 2 are good for 
detecting the scale-dependence of bias related to non-Gaussian processes. 

The statistics of the SFCs Wj^i are different from the WFCs Wj^i. As mentioned in §2.3, the 
scaling functions </>j,;(x) are window-like functions. SFCs are similar to the statistics of count in 
cell (CIC). The DWT scaling functions 0j,/(x) correspond to the CIC window functions, and the 
SFCs of the DWT correspond to the counts of CIC. Because 4)j^i{x) and the window functions of 
CIC are not orthogonal with respect to j (or scale.), the SFCs, u/^y, are given by a superposition 
of WFCs, Wji^v and / < j [see, Eqs.(8), (9) and (10)]. Namely, for any order of n, statistics based 
on SFCs are always sensitive to the location, or the phases of perturbations on different scales, i.e. 

in) 

sensitive to non-Gaussianity. Thus, regardless whether n is larger than 2, hj are always useful to 
detect the scale-dependence of bias related to non-Gaussian processes. 



3. Samples of galaxies 

Neither 6^"^ nor h^^'^ are directly testable because we don't know the distribution of dark 
matter bmix)- But, bias scale-dependence can be revealed from the difference between the scale 
behaviors of the galaxies with different internal properties. It has been very well known that the 
relative abundances of galaxies with different morphology are environment dependent (Dressier 
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1980), the clustering strength of galaxies is luminosity dependent (Xia, Deng &i Zhou 1987), and 

the correlation amplitudes of optically selected galaxies are different from that of infrared selected 
(Saunders, Rowan- Robinson & Lawrence 1992, Strauss et al. 1992). Therefore, it is most likely 
that galaxies with different morphology possess different scale-dependence of bias. Thus, the bias 
scale-dependence might be detected by comparing the j-spectra of 6" and bj for galaxies with 
different morphology. 

Let us consider two types of galaxies, I and II, both of which obey linear model (1) with 
different parameter b. All statistics of Eqs.(15) and (20) are available by replacing (5,m) by 
(7, II) . Namely, the scale dependence of bias can be revealed by j-spectrum of bj^^ and bj^^ 
defined as 

Any non-flatness of these j-spectra indicates that the bias is scale-dependent for at least one type 
of the considered galaxies. 

We analyzed the samples of galaxies listed in the APM bright galaxies catalog (Loveday, 
1996), which gives positions, magnitudes and morphological types of 14,681 galaxies brighter 
than 16"*. 44 over a 4,180 deg° area in 180 Schmidt survey fields of south sky. The completeness 
is about 96.3 per cent with a standard deviation 1.9 per cent inferred from carefully checked 12 
fields. Therefore, it is large and uniform enough for a 2-D DWT analysis. 

We choose the early type or elliptical and lenticular (EL) galaxies as /, and late type or 
spiral (SP) galaxies as II. All the elliptical and lenticular galaxies are compiled into one sample 
containing 4,439 ELs. The SP sample contains 8,217 SP galaxies. 

In order to do a 2-D DWT analysis, we chose three fields SI, S2 and S3 from the entire survey 
area. The three fields are selected to be square on a equal area projection of the sky. The equal 
area projection keeps the surface number density of the galaxies in the plane to be the same as on 
the sky. The sizes of SI, S2 and S3 are taken to be as large as possible covering the whole area 
of the survey. The whole survey and the three fields SI, S2 and S3 are plotted in Fig. 1. Each 
field has angular size of about 37° x 37°. SI, S2 and S3 contain 1,095, 1,039, and 1,055 ELs, and 
2,186, 2,188, and 2,092 SPs, respectively. The galaxies in the regions SI and S3 are completely 
independent. S2 has some overlaps with SI and S3. 

We divide each square region into 2^^^ x 2^^^ (1024^) cells labelled by 1 = (hjh), and Zi, I2 
being an integer from to 1023. Using the estimation of the mean depth of the sample given by 
the luminosity function from the Stromlo-APM redshift survey (Loveday et al. 1992), one find 
that the cell size is about 65 h~^Kpc, which is fine enough to detect statistical features on scales 
larger than 1 h~^Mpc. 
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The distribution of the galaxies in the sample can be described by a 2-D density field p{x.) 
or contrast (5(x), where x = (xi,X2)- In doing 2-D DWT analysis of /o(x) or 5(x), the wavelet 
functions and the scaling functions are constructed from direct product of the 1-D bases, i.e. 

= i^ji,h{^i)tpji,l2i^2), and ^/)j,i(x) = (I)ji,hi^i)(l)ji,i2{x2), where j = (ji,j2) and I = {li,l2)- 
For scale j, the angular scale is 2^'(^°~-') in unit of a cell. In this paper, we are interested in scales 
of j = 8,7,6,5,4, and 3, corresponding to spatial scales of about 0.26, 0.52, 1.04, 2.08, 4.16 and 
8.32 h~^ Mpc, respectively. 

In doing statistical analysis of the APM-BGC, the effect of the 1,456 holes drilled around big 
bright objects should be taken into account. For instance, random samples are generated in the 
area with the same drilled holes as that in the original survey. We also removed a few galaxies 
placed in the holes drilled. 

A common problem of statistics for discrete sample is sampling error, which is, in particular, 
serious for statistics of non-Gaussianity. Even if the original matter field is Gaussian, the sampled 
data must be non-Gaussian on scales for which the mean number in one cell is small. This 
is the non-Gaussianity of shot noise. Any non-Gaussian behavior of the density fields will be 
contaminated by the shot noise. Because of bias might be related to non-Gaussian process, the 
detected result is inevitably contaminated by sampling. Fortunately, the DWT spectrum method 
is found to be eff'cctive for suppressing the contamination of shot noise. It has been shown that 
the non-Gaussianity of shot noise is significant only on scales comparable with the mean distance 
of nearest neighbors of objects. Because the wavelet bases are localized, the central limit theorem 
guarantees the non-Gaussianity of shot noise rapidly and monotonously approaching zero on larger 
scales (Greiner, Lipa & Carruthers 1995, Fang &; Pando 1997). Therefore, besides small scales, 
the behavior of DWT j-spectrum will not be affected by the shot noise of sampling. 



4. Detection of j-spectrum 
4.1. Two point correlation function 

~(n) 

Before calculating bj , we consider the two-point angular correlation functions of the ELs 
and SPs of the APM-BGC. Because the two-point angular correlation function analysis of 
Stromlo-APM redshift survey has been done by Loveday et al. (1995). One can show the reliability 
of the DWT method by comparing our result with the earlier one. 

Fig. 2 shows the angular correlation functions, wel{0) and wsp{0), which are obtained from 
whole sample of the ELs and SPs. The dotted lines give the 1 a error calculated from 20 random 
samples which are produced by randomizing the positions of the galaxies in the same field with 
the same drilled holes as the samples. 

As usual, these correlation functions show a power law, and the amplitude of the correlation 
for the ELs is higher than that of the SPs. This is the well-known segregation of morphology: the 
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clustering of EL galaxies is stronger than that of SP (Davis & GcUcr 1976.) Actually, this is an 
evidence of bias. However, Fig. 2 shows that the ratio wel{&)/wsp{0) is approximately constant 
in the angular range of —0.4 < \og9 < 0.6. Namely, the segregation can be explained by linear bias 
if hi and bn are taken to be different constant. In other words, in terms of two-point correlation 
functions, no scale-dependent bias is needed, or at most a very weak scale-dependence. With the 
DWT analysis, we are able to measure this scale-dependence more quantitatively. 

This result of Fig. 2 will not be changed if the error is estimated by bootstrap samples. This 
is expected, because the number of galaxy pairs is large, the bootstrap-resampling error is only 
larger than Poissonian error by a factor of 1.7 (Mo, Jing & Borner 1992.) 



4.2. Statistics sf^ and sf^ 

Power spectrum is the Fourier counterpart of two-point correlation function. As mentioned 
in §2.4, the second order statistic {w^jj) is equivalent to the measure of Fourier power spectrum 

(Pando & Fang 1998). Therefore, one can expected that the second order statistics ^f^) should 
show the same results as that given by two-point correlation function. 

The 2-D extension of the n = 2 statistic {Vf) of Eq.(23) is 



r(2) ^ LZ^;i=o ^^2=0 "^j,iJ^ /2b^ 
[Z^/i=o l^i2=o 

In order to estimate the error of bj_^ and to avoid the effect of the drilled holes and boundary, we 

w 

3,3 



~(2) 

generated 500 randomized samples of the EL and SP, and calculated the average b, ■ and variance 



~ (2) 

from these random samples. Fig. 3 plots defined as 



5f = O^. (26) 

3,3 



~ (2) 

The variance from the random samples is also plotted in Fig. 3. The j-spectrum of B^- , indeed, 

~ (2) 

shows flat within l-cr, only field S2 has a little higher B^ on large scale j = 3. This result is in 
good agreement with the detection of two-point correlation functions (Loveday et al. 1995.) This 
shows that the DWT power spectrum estimator is reliable (Pando k, Fang 1998.) 

It has been emphasized in §2.4 that the flatness of the j-spectrum of the second order 

~ (2) 

statistics B)j doesn't imply that other j-spectra will also be flat, especially in the case that 
the bias scale-dependence is caused by non-Gaussian process. Therefore, non-Gaussian sensitive 

~ (3) 

statistics are necessary. Fig.4 shows the result of statistics B - which is defined by 

1(3) 

5(3) ^ (27) 
^ 1(3) ' ^ ' 

3,3 
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where 

E2J1-1 Y^2J2-1 I ~ |3i 
= <i=0 2^^2=0 l^j.ll \^ (28) 

31,32 r^2-'l-l Y^2^2-l I ~ 131 ■ ^ ' 

Comparing Fig.4 with Fig.3, it is clear that the j-spectra of show more deviation from a 
fiat spectrum. Therefore, the bias scale-dependence might really originate from non-Gaussian 
processes. 



4.3. Statistics Aj 

A simplest phase-sensitive statistic is the reconstructed distribution f>'{x) [Eq.(9)], which is 
given by SFCs. The 2-D extension of Eq.(9) is 

2n_i2-'2-i 

= E E «^j,i<^j,iW- (29) 

j = (ji,j2) means a reconstruction of the density field on scale ji in the dimension xi, and j2 in 
X2- Therefore, when j = (ji,0), /f>'(x) is a projection on axis x\, and when j = (0,^2), a projection 
on axis X2- In the case of j\ = j2 = j-, the field is smoothed on the scale j in both directions xi 
and X2, 

A DWT reconstruction on the three fields SI, S2 and S3 was performed. Since the SFCs, w^^i, 
are proportional to the density at position 1, the mean density of the reconstructed density field is 
proportional to 

I I 2n-l 2J'2-i 

m = —— y y w;i. (30) 

J 231 2^2 J'' ^ ' 

«1=0 «2=0 

The cells with Wj^\ > Wj are dense regions (clumps), and cells with w^^i < wj are underdense regions 
(voids). Actually, for all j, {wj)i and {wj)n are always proportional to the mean number density 
of galaxies EL and SP, respectively, and therefore, the ratio (wj) i / (wj) j j doesn't depend on j. 

Fig. 5 shows reconstructed distributions of SI on scales ji = j2 = j = 6,5,4,3. in which the 
bold contours denote the area with Wj^i-^-j^i^ — wj>2,3, 4...aj from outside to inside successivefully, 
and the 1 a line is calculated from 500 random samples for each field. Fig.5 also shows 
the underdense regions by light contours corresponding to wj^i^-j^i^/w^ < 1/2,1/8,1/32,1/128 
successivefully from out to inside. 

A simplified representation of these fields is shown in Fig.6, in which the dark area represents 
the overdense regions, Wj^i^^-j^i^ > wj, and blank the underdense regions, wj^i^-j^i^ < wj. The ratio 
between the clustering strengths of the EL and SP on scale j can be measured from Fig.6 by 

^ _ (dark area)g/, 
(dark area)5P 
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Similar to the statistic Ir- ' , to estimate the error of aj , we generated 500 randomized samples of 
the ELs and SPs, and calculated the average aJ of these random samples. Fig. 7 plots the result 
of statistics Aj defined as 

A, = X (32) 

Uj 

The one a lines of both randomized and bootstrap re-sampling are also shown in Fig.7. The 
two error estimates gave about the same confidence level. The j-spectrum of Aj is found to be 
significantly non-flat. 



4.4. Statistics Bj 



Since Aj is based on statistic Wj_\, one can expected that the statistic of b^^^ of Eq.(21) will 
give the same conclusion. The 2-D extension of Eq.(21) is 



E 231-1 ^2^2-1 1 11 



where we use absolute value \w^^\ — because for odd n, {{w-^^\ — w-^)'^) has larger relative errors. 
We calculated B^p defined as 

Bf = (34) 

where 6^^ is the average of j over 500 randomized samples of the ELs and SPs. The results of 

B^P for the three fields are given in Fig. 8. The three fields SI, S2 and S3 show the same feature 
of the j dependence. The mean SFCs and their variance from the randomized samples are also 
plotted in Fig.8. Like the statistics Aj, the spectrum of b'^P is non-flat with significance larger 
than 2 a. 

Since 2-D sample of galaxies is a z (redshift) projection of their 3-D distribution. To study 
the influence of sample depth on our statistics, we calculated B^ for samples of ELs and SPs 
with different limit magnitude m. Fig.9 plots the result of -B]^^^ which is the same as B^ but for 
sample with limit magnitude m = 16.20. In the case, the numbers of EL and SP galaxies are about 
20-25% less than that of m = 16.44. The j-spectra in Fig.9 have completely the same features as 
Fig.8. Therefore, conclusions drawn from Fig. 8 are insensitive to the depth of the samples. 

In Figs.8 and 9, the j-dependence of bias is prominent on scales larger than 1 h""^ Mpc, i.e. 
larger than the mean distance of nearest neighbors of galaxies, and therefore, the effect of shot 
noise is negligible. 
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5. Discussions and conclusions 

We showed that, instead of one parameter b, Eq.(l) introduces a series of parameters 6^"^ and 

bj to describe the bias scale-dependence. The statistics of 5^- with different n are independent 
from each others. The DWT analysis provides a simple and effective tool of systematically 
detecting the scale-dependence of bias parameters on various orders. This method is effective to 
be employed for analyzing galaxy samples which show morphology- and/or luminosity-segregation. 

~(n) (n) 

With this method we detected the j-spectra of bj and bj for the distributions of EL 
and SP galaxies of the APM-BGC samples. The general results indicate that for second order 
statistics, i.e. two-point correlation function and power spectrum, the bias is approximately 
scale-independent, but not so for higher order or phase-sensitive (non-Gaussian) statistics. The 
result is consistent with the following fact: most evidences for weak scale dependence of bias are 
from statistics of the two-point correlation functions and power spectrum (Kauffmann, Nusser & 
Steinmetz 1997), while the evidence for scale-dependence is from phase-sensitive statistics (Sigad 
et al. 1998). Therefore, the scale-dependence of galaxy bias may have the same origin as the 
non-Gaussianity of galaxy distribution. Linear evolution cannot cause non-Gaussianity of mass 
distributions if the initial perturbations are Gaussian. The scale-dependence of bias of galaxy 
distribution is most unlikely due to the non-linear evolution of gravitational clustering, and the 
non-local relationship between galaxy formation and their environment. 

With this in mind, the information of bias scale-dependence is worth for developing models 
of galaxy formation. Indeed, despite the current detection of bias scale dependence is still very 
preliminary, the result is already able to set useful constraint on models of hierarchical clustering. 
In these models, galaxy correlations are generally assumed to be described by the hierarchical 
relation ^„ = QnC2^^ where ^„ is the n-th order correlation function, and Qn are constants 
(White 1979). If these hierarchical relations hold exactly, the second order (two-point) correlation 
function plus all constants Qn (which may be different for different types of galaxies) completely 

~ (n) 

characterize the clustering of galaxies, including their higher order correlations. This is, all Bj 
and Bj^^ can be represented by second order correlation function plus all (scale-independent) 
constants Qn- Hence, if bias is scale-independent on second order, it will be scale-independent on 

all orders. Therefore, the non-flatness of j-spectra of Bj^^ and Bj^^ for sample APM-BGC implies 
that the hierarchical relations may not hold exactly, or the coefficients, Qn, are scale-dependent. 
Similar conclusion has also been drawn from the detection of the scale-scale correlations of the 
Lya forests of QSO absorption spectrum (Pando et al. 1998.) Thus, the detection of sj"^ and 

(n) 

B^ , joining with other higher order statistics, is effective to reveal the details of the hierarchical 
clustering scenario. 

We thank Drs. J. Pando and Y.P. Jing for many helpful comments. ZGD and XXY were 
supported by the National Science Foundation of China. 
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Fig. 1. — The equal project of the APM bright galaxy distribution and three regions SI, S2 and 
S3, in which the wavelet analysis has been done. The top and bottom panels are for EL and SP 
galaxies respectively. 
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Fig. 2. — Two point angular correlation functions of EL (solid lines) and SP (dashed lines) galaxies 
in APM bright galaxy catalogue. 
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Fig. 3 



~ (2) ~ (2) 

Fig. 3. — Bj vs. j for three fields SI, S2 and S3. The mean of and 1 a error given by 
randomized samples are plotted by dot-dashed and dashed lines, respectively. 
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Fig. 4 



~ (3) ~ (3) 

Fig. 4. — Bj vs. j for three fields SI, S2 and S3. The mean of and 1 a error given by 
randomized samples are plotted by dot-dashed and dashed lines, respectively. 
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Fig. 5 

Fig. 5. — Reconstructed distributions of SI on scales ji = j2 = 6, 5,4, 3 from top to bottom. Left 
side results are for EL galaxies and right for SPs. The bold contours are for Wj^i^-j^i^ —wj > 2, 3, 4...aj 
from outside to inside successivefuUy, and the light contours denote the area with wj^ij^.j^i^/w^ < 
1/2,1/8,1/32,1/128 from outside to inside successivefuUy. 
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Fig. 6 

Fig. 6. — A simplified representation of Fig. 5. Dark area denote the overdense regions, 
Wj^l^-j^l^ > Wj, and the blank area the underdense regions, wj^i^-j^i^ < wj, where ji = j2 = j- 
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Fig. 7 

Fig. 7. — Aj vs. j for three fields (solid lines). It shows the mean of Aj (thick solid lines), and 1 a 
errors given by randomized samples (dotted lines) and bootstrap-resampling (thin solided lines.) 
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Fig. 8 

Fig. 8.— Bf' vs. j for three fields (solid lines.) The mean of sf' and 1 a error given by randomized 
samples are plotted by dot-dashed and dashed lines, respectively. 




Fig. 9 



Fig. 9. — B^- vs. j for three fields (solid lines.) The mean of B^- and 1 aj error given by 
randomized samples are plotted by dot-dashed and dashed lines, respectively. 



